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A New Luminous Blue Variable in M31^ 
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ABSTRACT 


We report the hfth conhrmed Luminous Blue Variable/S Doradus variable 
in M31. In 2006, J004526.62+415006.3 had the spectrum of hot Fe II emissio n 
line star with strong P Cygni prohles in the Balmer lines flMassev et al.l 120071) . 
In 2010, its absorption line spectrum resembled an early A-type supergiant with 
H and Fe II emission lines with strong P Cygni prohles, and in 2013 the spec¬ 
trum had fully transitioned to an F-type supergiant due to the formation of the 
optically thick, cool wind which characterizes LB Vs at maximum light. The pho¬ 
tometric record supports the LBV/S Dor nature of the variability. Its bolometric 
luminosity ~ -9.65 mag places it on the HR Diagram near the known LBVs, AE 
And, Var C in M33 and S Dor. 


Subject headings: stars 


1. Introduction 


The class of stars known in the literature as Luminous Blue Variables(LBVs) are an 
important phase in the hnal evolution of evolved massive stars. They may represent a short 
high mass losing stage in which the star sheds much of its outer envelope prior to becoming a 
Wolf-Rayet star, or as some authors have suggested, they may be the immediate progenitors 
of some types of Type II supernovae. A classical LBV or S Dor variable in quiescence 
(minimum light) resembles a moderately evolved hot star usually with the spectru m of a B- 
type supergiant or an Of/late WN star (iHnmphreys fc Davidsorjll994t IVinki 120121 ). During 
an LBV eruption, the mass loss rate increases, the wind becomes opaque and cool, and 
its absorption spectrum resembles an A to F-type supergiant. Since this is a shift in the 
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bolometric correction, the star brightens in the visual and appears to move to the right, 
to lower temperatures on the HR Diagram. This is the LBV’s optically thick wind stage 
or maximum light. The origin of the instability that leads to this dramatic change is not 
understood, but most suggestions involve an instability close to the Eddington Limit. The 
mass loss event can last for several years or even decades, and afterwards the star returns to 
it previous state. 


Many stars are called LBVs or candidate LBVs, but few are confirmed probably due 
to the infrequency of the LBV “eruption” or maximum light stage. For e xample, there are 


only f our conhrmed classical or normal LBV/S Dor variables in M31 f.see iHnmphrevs et al. 


(120141) ). Thus the discovery of a new LBV in its eruption state is important for understanding 
their relation to the larger massive star population, the frequency of the LBV eruption, and 
the duration of this unstable phase. 

J004526.62+415006. 30 was described by Massev et ah f 2007|) as a “hot LBV candidate” 
Its normalized spectrum from 2006 (Fig 12 in iMassev et al.l (120071 1) showed strong H emission 
lines with P C ygni profiles and weak F e II emission lines. In a survey of luminous stars in 
M31 and M33. iHnmphrevs et ahl (120141) noted t hat its 2010 spectrum clo sely resembled that 
of the warm hypergiant J004444.52+412804.0 (IHnmphrevs et ahl 120131) with prominent P 
Cygni prohles in the multiplet 42 Fe H emission lines, strong hydrogen emission with broad 
wings and P Cygni prohles, and the absorption line spectrurn of an early A-type supergiant. 
A spectrum from 2011 published bv ISholukhova et al.l (120151) clearly showed that the Fe H 
lines had weakened. They also noted the change from 2006 and suggested that it is an LBV. 
A second star J004051.59+403303.0 was identihed as an LBV candidate. 


In this paper we discuss recent spectra and photometry of J004526.62+415006.3 that 
conhrm it is an LBV currently in an optically thick wind state or maximum light. 


2. New Observations 


In addition to the spectrum described in Humphreys et ahl (l2014l) . it was observed 
again on 25 and 26 S e ptern ber 2013 with the Hectospec Multi-Object Spectrograph (MOS) 
( Fabricant et al 1998 . 2005 ) on the 6.5-m MMT on Mt. Hopkins. The Hectospeca has a 1° 
FOV and uses 300 hbers each with a core diameter of 250//m subtending 1'.'5 on the sky. We 
used the 600 1/mm grating with the 4800A tilt yielding ps 2500A coverage with 0.54 A/pixel 


^M31-004526.62 in Humphreys et al. ( 2014h and Sholukhova et al. (2015) 
^ http://WWW. cfa. harvard. edu/mmti/hectospec. html 
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resolution and R of ~ 2000. The same grating with a tilt of 68 OOA was used for the red 
spectra with 2500A coverage, 0.54A/pixel resolution and R of ~ 3600. The spectra were 
reduced using an exportable version of the CfA/SAO SPECROAD package for Hectospec 
data§|. The spectra were all bias subtracted, flathelded and wavelength calibrated. 

Recent broadband CCD photometry was obtained by J. C. Martin with the 20-inch 
telescope at the Barber Observatory in 2013 and 2014, and by Sholukhova et ah (2014) 
in 2011. We also measured additional magnitudes from /f^T/ACS/HRC images observed 
August 2010 (HST-GO-12056) and January 2012 (HST-GO-12106). The magnitudes are 
zero pointed to Vega and comparable to Johnson B and Cousins I. 

All of the spectroscopic and photometric observations discussed here are summarized in 
Table 1. 


3. The Spectrum from 2006 to 2013 and Light Curve 


Phil Massey kindly provided his spectrum from 2006 flMassev et al.l 120071) reproduced 


here in Figure 1. iMassev et al.l ( 20071) described J004526.62-1-415006.3 as a “hot LBV can¬ 
didate” . The spectrum is dominated by strong hydrogen emission with broad wings and P 
Cygni prohles and Fe II and [Fe II] emission. Some weak absorption lines of He I A4026, O 
II A4070-4076 and S IV A 4116 are visible, suggestive of an early B-type star. 


The blue and red spectra from 2010 and 2013 are shown together in Figures 2a and 2b, 
respectively. They not only illustrate the dramatic change from the hot emission line star in 
2006 but also the transition to a cooler temperature from 2010 to 2013 and t he formation 
of the optically thick cool wind in the 2013 spectrum. iHnmphrevs et al.l (120141) classihed its 
2010 spectrum as an early A-type supergiant (A2ela). With its strong P Cygni prohles in 
the multiplet 42 Fe II e mission lines they noted its close resemblance to the warm hypergiant 
J004444.52-1-412804.0 ( Humphreys et ah 2013). 


The 2013 spectra show the continued evolution of the wind to cooler temperatures. 
Metallic absorption lines have replaced the Fe H emission spectrum, strong Ca H H and K 
lines are conspicuous, and the P Cygni prohles are gone except for Ha. The blue absorption 
line spectrum of the optically thick wind is consistent with an early F-type supergiant 
F2). In the red spectrum, the luminosity sensitive O I A7774 triplet and the N I (multiplet 
3) absorption lines are much stronger. The diherences between the 2010 and 2013 spectra 


^External SPECROAD was developed by Juan Cabanela for use on Linux or MacOS X systems outside 
of CfA. It is available online at: http://iparrizar.mnstate.edu. 
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Wavelength 


Fig. 1.— The Fe II and hydrogen emission line spectrum from 2006 flMassev et ahll2007() 


illustrate the spectroscopic transition and the formation of an optically thick wind typical 
of LBV/S Dor variables. 


Like other LB Vs I Hnmphrevs et al.l 120141: lOksala et al.! 120131: iKraus et al.l 120141) . its 
spectral energy distribution has a small near-infra red excess due to free-free emission, but 
it was not detected in the IRAC survey of M31 ( Mould et ah 2008) or the WISE Source 
Catalog. So it apparently do e s not have a near- or mid-infrared excess due to dust (see 
Figure 11 in iHnmpIirevs et ahl (2014)). However, the relative strengths of the [N II] A5755A 
and A6584A emission lines compared to the [S II] lines in the 2010 spectrum indicate the 
presenc e of a relatively d ense circumstellar nebula, much denser than expected for an H II 
region (jWeis et aLll2014l ). The outflow velocities from the stellar wind measured from the 
absorption minima in the P Cygni prohles in these spectra are summarized in Table 3. LBVs 
not only have slow wind velocities during their maximum light, dense wind state, but during 
their quiescent hot state their winds are als o slower than those me asured the same way in 
normal supergiants of similar spectral type flHnmphreys et ahl 120141 1. J004526.62-1-415006.3 
shows the same pattern, although the outflow measured in the 2010 spectrum is interestingly 
higher when the star was in transition to the dense wind state. 


The light curve is shown in Figure 3. Unfortunately, we do not have either spectra 
or photometry between 2006 and 2010 to determine the onset of the LBV/S Dor eruption. 
However, our spectra demonstrate that the star was in transition to its cool dense wind 
state in 2010, so the outburst probably began in 2010 or shortly before that. Based on the 
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Fig. 2a.— The blue spectra from 2010 (blue) and 2013 (red). Note the development of the 
absorption line spectrum replacing the Fe II emission lines and weakening of the P Cygni 
prohles. 


photometry from 1999-2003 to 2012-2014, AV is surprisingly, at most, only 0.9 mag which is 
too small for the implied change in apparent temperature between the 2006 and 2010-2013 
spectra. Thus we suspect that the earlier photometry is not representative of its minimum 
light or quiescent state, and it may have been in an elevated state at that time. Previous 
photometry also show s that it has had other eruptions; it was measured at V = 16.3 in 1992 
f Magnier et ah 1992). 


4. Relation to Other LB Vs 


Adopting the interstellar extinction, = 1.5 mag from iHumphrevs et al.l (120141) . mea¬ 
sured from three different indicators, its visual luminosity at rn aximum light is -9.65 
with a distance modulus of 24.4 mag for M31 flRiess et al.ll2012l) . At maximum light, the 
bolometric correction will be quite small, corresponding to the early F-type spectrum of the 
dense wind. is therefore also ~ -9.65 mag. The corresponding parameters for the star 
in its quiescent or hot state are less well determined, but its spectrum from 2006 suggests an 
early B type supergiant with a surface temperature on the order of 20,000 K or hotter. On 
the HR Diagram, J004526.62-1-415006.3 will thus he very close to AE And in M31, Var C in 
M33, and the well studied progenitor of the class S Dor, see Figure 12 in [Humphreys et ah 
(l2014l) . 
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Fig. 2b.— Red spectra from 2010 (blue) and 2013 (red). The luminosity sensitive O I triplet 
at A7774A in intermediate temperature supergiants is much stronger in 2013. 

Given the lack of photometric data between 2000 and 2010, any additional data from 
that period or earlier would be appreciated to better define the properties of the progeni¬ 
tor and the onset of the current eruption, as well as future photometric and spectroscopic 
monitoring. 

We are especially grateful to Phil Massey for sharing his spectrum of J004526.62-1-415006.3 
from 2006 and to Olga Sholukova and her collaborators for a preprint of their paper. We 
also think the Hectospec snpport operators Perry Berlind and Mike Calkins. Research by R. 
Hnmphreys on massive stars is snpported by the National Science Fonndation AST-1109394. 

J. C . Martin’s collaborative work on Inminous variables is snpported by the National Science 
Foundation grant AST-1108890. 

Facilities: MMT/Hectospec, HST/ACS/HRC 
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Fig. 3.— The light curve based on the the data in Table 2. The symbols are B, open blue 
squares; V, hlled green circles, R, red crosses and I, open red triangles. A gap is show n 
between 2000 and next available data in 2010. 
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Table 1. Journal of Observations Used in this Paper 


UT Date 

Instrument 

Exp. Time 

Wavelength or Filters 

Resolution 

Reference 

Spectroscopy 

19-20 Sep. 2006 

WIYN/Hydra 

180m 

3970-5030A 

i.sA 

Massev et al. ('20071 

21-22 Sep 2006 

WIYN/Hydra 

90m 

4550-7400A 

3.4A 

Massev et al. ('20071 

10 Oct. 2010 

MMT/Hectospec 

120m 

3600-6050A 

1.9-2.2A 

Humphreys et al. 2013,2014) 

09 Oct. 2010 

MMT/Hectospec 

90m 

5500-8000A 

1.9-2.2A 

Humphreys et al. 2013,2014) 

25 Sep. 2013 

MMT/Hectospec 

120m 

3600-6050A 

1.9-2.2A 

This paper 

26 Sep. 2013 

MMT/Hectospec 

90m 

5500-8000A 

1.9-2.2A 

This paper 

Photometry 

1999 - 2003 

WFC/INT 


BV 


Vilardell et al. 120061 

03 -06 Oct. 2000 

Mosaic CCD/4m 


UBVRI 


Massev et al. 120061 

13 Aug. 2010 

HST/ACS 


BI 


This paper 

Oct.-Nov. 2011 

BTA 6-m. 


BVR 


Sholukhova et al. 120151 

11 Jan. 2012 

HST/ACS 


BI 


This paper 

02 Sep. 2013 

CCD/Barber Obs. 


VR 


This paper 

28 Oct. 2013 

CCD/Barber Obs. 


VR 


This paper 

21 Oct. 2014 

CCD/Barber Obs. 


BVR 


This paper 

21 Nov 2014 

CCD/Barber Obs. 


BV 


This paper 

15 Jan 2015 

CCD/Barber Obs. 


BVR 


This paper 
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Table 2. Multi-color Photometry 2000 - 2015 


Date 


1999 - 2003 
03 - 06 Oct. 2000 
13 Aug. 2010 
Oct.-Nov. 2011 
11 Jan. 2012 
02 Sep. 2013 
28 Oct. 2013 
21 Oct. 2014 
21 Nov. 2014 
15 Jan. 2015 


B mag 

V mag 

17.11 

17.08 

17.66 

17.16 

17.31 


16.82 

16.37 

16.81 

16.23 


16.27 

16.66 

16.46 

16.68 

16.45 

16.76 

16.51 


R mag 

I mag 

17.48 

17.35 


16.61 

16.08 

15.93 

15.93 


15.99 


16.17 


16.27 



Table 3. Outflow Velocities in J004526.62-1-415006.3 


Date 

Velocity km s ^ 

Lines 

Sep. 2006 

-101 ± 3 

3 Hydrogen lines 

Oct.2010 

-149 ± 6 

4 Hydrogen lines 


-138 ± 3 

3 Fe H lines 

Sep.2013 

-115 

Ha only 













